NONATTACKING QUEENS IN A RECTANGULAR STRIP 



SETH CHAIKEN, CHRISTOPHER R.H. HANUSA, AND THOMAS ZASLAVSKY 

Abstract. The function that counts the number of ways to place nonattacking identical 
chess or fairy chess pieces in a rectangular strip of fixed height and variable width, as a 
function of the width, is a piecewise polynomial which is eventually a polynomial and whose 
f*^) ' behavior can be described in some detail. We deduce this by converting the problem to one 

£NJ , of counting lattice points outside an affinographic hyperplane arrangement, which Forge 

and Zaslavsky solved by means of weighted integral gain graphs. We extend their work 
. by developing both generating functions and a detailed analysis of deletion and contraction 

for weighted integral gain graphs. For chess pieces we find the asymptotic probability 
that a random configuration is nonattacking, and we obtain exact counts of nonattacking 
<N ■ configurations of small numbers of queens, bishops, knights, and nightriders. 
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^ ; 1. Nonattacking pieces 

H ■ 

A classic enumeration problem is that of counting the configurations of eight nonattacking 
queens on a chessboard, and more generally of n nonattacking queens on an n x n board. 
While thinking about how to represent this problem geometrically we realized that a similar 
problem, a generalization in which the height of the board is fixed but the board's width 
can vary, is amenable to treatment by a recent geometrical counting method of Forge and 
Zaslavsky [2]. 

Here is our variant question: Place m nonattacking chess queens in a rectangular strip of 
fixed height m and variable width n. In how many ways can this be done? We show that, 
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as n increases, the answer v m {n) becomes a polynomial in n. Indeed, the number of ways is 
a sum of more and more terms as n increases from 0, each term being zero until n reaches 
a threshold after which the term is a product of linear factors. There is one term, n m , with 
largest degree, and this term is present for all nonnegative n; thus one may say that the 
number of nonattacking placements is a piecewise polynomial of degree m. Furthermore, 
the signs of the terms are such that there is no cancellation between terms. All this applies 
to any chess piece with centrally symmetric moves, placed one to a row; that is, to all but 
the pawn. Indeed it applies to any fairy chess piece with suitable moves. A fairy chess piece 
is an imaginary chess piece with an arbitrary movement rule; many such pieces have been 
tried and some found exciting; we mention the queen with a knight's move added to its usual 
repertory, or the nightrider, which extends a knight's move indefinitely. 

Forge and Zaslavsky counted lattice points inside a hypercube of variable side length but 
not lying in any of a finite set of affinographic hyperplanes (hyperplanes determined by 
constancy of the difference of two coordinates), by converting the problem into one about 
weighted integral gain graphs (which will be defined later). Their principal result is that the 
counting function is a piecewise polynomial with predictable leading term, and is eventually 
a polynomial. This yields our basic theorem because, in the space of configurations of chess 
pieces, attack is an affinographic property. 

From knowing the moves of our piece we also get more detailed information about the 
counting function. For pieces like the queen, bishop, and nightrider, we find the exact value 
of m at which the function becomes a polynomial (Proposition [3]). When there is one piece 
in each row, the magnitude of the second coefficient of the counting polynomial has a simple 
combinatorial interpretation; it equals half the maximum total number of squares attacked 
nonhorizontally if the pieces are placed on a wide board and they do not attack any of the 
same squares (this is a special case of Proposition [6]). This second coefficient is important 
because it determines the asymptotic probability that a random placement of one piece in 
each row will be nonattacking (Corollary [7j). 

Besides this, we extend the general theory of weighted integral gain graphs to treat gener- 
ating functions and the details of the method of deletion and contraction; then we apply the 
latter to get exact formulas and generating functions for the number of nonattacking con- 
figurations with fixed small height for four chess and fairy chess pieces: the queen, bishop, 
knight, and nightrider. We do not find any general formula for the counting function; that 
seems impractically difficult, on the order of finding the chromatic number of a graph. 

2. Fundamentals 

Here is the exact setup in complete generality. We have a rectangular board with m rows 
and n columns, m being a fixed and n a variable positive integer. We think of the board as 
the set [n] x [m] C M 2 , where [n] '■— {1,2, . . . ,n}. We have a chess piece P — a fairy chess 
piece, which means its moves are arbitrary (subject to two rules to be stated); the set of 
moves is Ap C Z 2 and signifies that if P is located at a position (j, k) G Z 2 , then it attacks 
all the squares in (j,k) + Ap := {(j,k) + (h,i) : (h,i) G Ap} that lie in the board. (We 
say "moves to" as a synonym of "attacks" ; the actual moves, if different from the attacked 
squares, as with a pawn, are irrelevant to our problem.) We need not assume (0,0) G A P ; 
but we do have two requirements for the moves. First, only finitely many moves can have 
any one nonzero value of the vertical displacement. Second, the set of horizontal moves, 
{(h, 0) G Ap}, should be either finite or the infinite set Z x {0}. We shall also assume that 
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Ap is centrally symmetric, but this is just a technical convenience to ensure that if one piece 
(a copy of P) attacks another, the attack is reciprocal, and we can ensure it, if necessary, 
through replacing Ap by Ap U (—Ap). Besides this, for each row j we specify a number qj 
which is the exact number of pieces to occupy that row. (The case qj > 1 is possible only 
if P has a finite set of horizontal moves. Then the "move" (0, 0) is important. It prevents 
two pieces from occupying the same square, and if it is not in the move set, the number of 
pieces on a square is limited only by the number of pieces in the row.) 

Thus, we can describe the entire configuration of q := qi + q 2 + - ■ - + q m pieces by q variables 
Xp where 1 < j < m and 1 < r < qj and 

x r j = the column occupied by the rth piece in row j. 

We call the vector x := (a£) a labelled configuration. Since the pieces are identical but the 
setup treats them as distinguishable, two labelled configurations describe the same configu- 
ration of pieces if they have the same occupied positions, that is, the same multiset of values 
of the variables in each row. Thus, the number of nonattacking configurations, v m {n), is the 
number of labelled configurations, \ m (n), divided by the number of labellings: 

(1) u m (n) 



Treating a labelled configuration as a point in [n] q C M 9 , the rules of nonattack are simple. 
A configuration is nonattacking precisely when: 

(Nl) if 1 < i < j < m, then (x| — x\,j — i) A P for all r G [q^\ and s G [qj], and 
(N2) when r < s in [qj], then (x* — x r j7 0) ^ A P . 

Rule (Nl) is a finite list of requirements because j — i is bounded by m — 1. Write 
j — i = k 7^ 0. Suppose there are 4 moves with height k and they are (fiki, k), . . . , (/ifct fe , k). 
(Possibly t k = 0.) Then (Nl) says 

x s i+k ^xl + fi k i 

for / G [tk] and ie [m — k]. 

Rule (N2) is also finite. If the set of horizontal moves is finite, say ±(/i i, 0), . . . , ±(// t , 0), 
(N2) says 

x S j ^ x] ± [i i 

for / G [to]- If it is infinite, then it is Z x {0} so (N2) simply requires that all qj < 1. 

Theorem 1. The number of nonattacking configurations of pieces having qj pieces in row j 
for each j G [m] is a piecewise polynomial function ofn and is a polynomial ifn is sufficiently 
large. The piecewise polynomial has degree q = ^2qj and leading coefficient (qi^- ■ • • qrn-)^ 1 ■ 
A sufficient condition for n to be sufficiently large is 

n > (q — 1) ■ max{|/i| : (//, k) G A P , < k < m}, 

(but we take 1 < k < m if P has unbounded horizontal moves or if all qj < 1). 

Later on, in Corollary \5\ we get a more detailed description of this counting function. 

For a piece with bounded horizontal moves, the total number of nonattacking configura- 
tions of q pieces is a piecewise polynomial given by summing the values given by Theorem [T] 
for all weak compositions q% + ■ — I- q m of q. Thus: 
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Corollary 2. The number of nonattacking configurations of q pieces with a finite number of 
horizontal moves is, for large n, a polynomial which has degree q and leading coefficient 



The simplest case is that in which every row has one piece, i.e., all qi — 1 and q = m. 
Then the factor gilq^! • • • (ZrJ = 1, v m (n) = A m (n), and the piecewise polynomial is monic. 
We solved this case for four different pieces, the queen, bishop, knight, and nightrider (a fairy 
chess piece with indefinitely extended knight's move) for small numbers of rows; see the last 
section for the results. The calculations quickly become too complicated for hand solution. 
Hanusa has prepared a pair of computer programs that manipulate weighted integral gain 
graphs to produce exact answers and generating functions; they are available on the Web [I] . 
The first phase, in Java, carries out operations on weighted integral gain graphs. The second 
phase, in Maple, uses John Stembridge's symmetric function package SF [6] to produce 
generating functions. The calculations uncovered an unexpected fact: in the polynomial 
that gives the number of configurations for large n, the second coefficient has a simple 
combinatorial meaning. See Proposition |6j 

3. Gain graphs, affinographic hyperplanes, and the integral chromatic 



An affinographic hyperplane in M. q is a hyperplane whose equation has the form Xj = Xi + c; 
it is called integral when c is an integer. An arrangement of hyperplanes is a finite set, 
considered in conjunction with the way it decomposes R q into <i-dimensional pieces (called 
regions of the arrangement). 

An integral gain graph $ consists of an underlying graph ||$||, with vertex set V = 
{v i, . . . , v q } and edge set E, and a gain function (p, which assigns an integer to each oriented 
edge e, with the rule that ^(e -1 ) = — (p(e) where e _1 denotes the same edge e with the 
opposite orientation. There are two kinds of edge: a link has two different endpoints, while 
a loop has both ends at the same vertex. The gain of a path P, ip(P), is the sum of the gains 
of the edges, which must be oriented in a consistent direction along the path. For brevity 
we write fiVjVji for an edge with endpoints Vj,Vji and gain /i in the direction from Vj to ty. 

In [2] integral gain graphs were used to treat integral affinographic hyperplane arrange- 
ments. To each edge e, with endpoints Wj and Vj, there is a corresponding affinographic 
hyperplane in W q whose equation is Xj = Xj + (p(e) if we calculate the gain with e oriented 
from Vi to Vj. This hyperplane does not depend on the orientation of e because, taking 
the opposite orientation e _1 from Vj to Vi, one gets the same equation. Thus every integral 
gain graph has a corresponding affinographic hyperplane arrangement, and conversely every 
integral affinographic hyperplane arrangement has a corresponding integral gain graph. 

We define xf ( n )> the integral chromatic function of $, as the number of mappings x : 
V — > [n] such that Xf ^ Xj + /i for each edge /iVjVf. (We may think of a function x as an 
integer lattice point in the hypercube [n] m that is not in any of the affinographic hyperplanes 
that correspond to the edges of $.) Let 



the largest gain of any path. According to [21 Corollary 3.2 and the following paragraph], 
x|(w) is a polynomial on the domain n > no($) but not on n > no($) — 1. 




qi >o 



FUNCTION 



rio(<fr) := ma.x{(p(P) : P is a path in $}, 
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To calculate by deleting and contracting edges, as we do in Section [51 we need 

an added feature: vertex weights. A weighted integral gain graph (<£>, h) (which we shall 
herein call a "gain graph" for short) is an integral gain graph $ which also has an integer 
weight hi assigned to each vertex t>j. The integral gain graphs introduced earlier, and any 
integral gain graphs without explicit weights, implicitly have weight at every vertex. Given 
a weighted integral gain graph, x% m ( n ) denotes the number of functions x : V — > Z such 
that Xj 7^ Xi + ^p(f) for each edge / from Vi to Vj, and hi < Xj < n for every vertex. We call 
xf§ ft ) the integral chromatic function of the gain graph. 

The deletion-contraction identity to be stated in Equation ([3]) later in the paper requires 
definitions of deletion and contraction of an edge e in a weighted integral gain graph. The 
definition of deletion is that e is removed from the graph with no other change; we write 
($, h)\e for the deleted gain graph. Contraction is more complicated. We write ($, h)/e for 
the contracted gain graph; but for computations a simpler notation is ($", h") := (<&,h)/e, 
in which the underlying graph is (V", E") and the gain function is ip". 

Before defining contraction we must define switching. Write V = {t>i, . . . , v q }. A function 
rj : V — > Z is called a switching function. Switching by rj changes the gain graph ($, h), with 
gain function if and weight function h, to ($ ,? , h v ) with the same underlying graph but with 
gain function <f v (f) : = <p(f) — ?7i + r/j, where / is an edge oriented from Vi to Vj, and with 
weights h\ := h k + r] k . 

Now fix an edge e with endpoints Wj and Vj and suppose it has nonnegative gain <p(e) 
when oriented from to Vj. One may think of contraction as a two-step process wherein 
one modifies Vi and its incident edges and then merges Vi into Vj. To do this we define a 
switching function t] by % := for k ^ i and rji := <f(e). To contract e, 

(1) first switch ($, h) by i], which changes the gain of e to and that of each other edge / 
to (f"(f) ■= f(f) — ^(e) if / is a link incident with Vi (we consider e, / to be oriented 
away from vi) and to f"(f) '■= f(f) otherwise; the weight of every vertex stays the 
same except for hi, which changes to K\ = hi + (p(e); 

(2) then coalesce the endpoints of e to form a new vertex v"j G V" and delete e, retaining 
the switched gains of all edges f ^ e and the weights of all vertices v u ^ Vi, Vj] and set 
the weight h"j of the contracted vertex v"j equal to max(^, hj) = max(/ij + <p(e), hj). 

(All this is from [2, Section 2], except that we have updated [2]'s rooted integral gain graphs 
to [5]'s more versatile weighted integral gain graphs.) Parallel edges with the same gain can 
be combined into one, because they represent the same constraint on coloring. 



4. Formulas for non-attacking configurations 

The particular integral gain graph $ we need for the chess problem has q vertices, one for 
each piece, labelled Vji for 1 < j < m and / G [qj] (representing piece (j, I), i.e., piece I in row 
j). There is an edge called iiVjiVyv for each I G [qj], I' G [qj'], each distinct j and j' = j + k 
in [m], and each /x such that (fi,k) G Ap (representing a possible attack (fM,k) from piece 
(j, I) to piece (f,l')). There is also an edge called iiVj\Vjy for each I ^ V in [qj] and each \x 
such that (/i, 0) G A P (i.e., for each potential attack [fi, 0) from piece I to piece I' in row j). 
The number /i is the gain of the edge; that is, <p(/iVjiVji>) = fi. These edges are directed but 
they always come in opposite pairs, fj,VjiVjqi and {—^VjiyVji, one should consider each pair 
as a single undirected edge. 
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Define hyperplanes in R 9 , with coordinates x r p by the equations 
(HI) x s i+k - xl = n H for k G [m- G [m - k],l G [t k ],r G [q t ],s G [ft +fc ], 
and, if is horizontally finite, 

(H2) - = ±Moi for j G [m], 1 < r < s < / G [t Q ]. 

The set "K of all these hyperplanes is the arrangement of affinographic hyperplanes that 
corresponds to the integral gain graph $ of the preceding paragraph. 

Proof of TheoremUl A nonattacking labelled configuration is a point x G [n] q \ (U JCj , and 
conversely. The main result of [2], Theorem 3.1, gives a formula for the number of such 
points, our \ m {n). The formula shows that X m (n) is a monic polynomial in n of degree equal 
to the dimension of the hypercube [n] q , if n is large. Moreover, [2, Corollary 3.2] shows how 
X m (n) behaves when n is not large; it is A m (n) = where 

ro d r 

(2) xl(n) = J2(- 1 ) 9 ' drc r l[max(0,n-n ri ), 

r=l i=l 

in which ro is a positive integer, the n r i are nonnegative integers, d\ — q and the first term 
is n q , the other d r < q, and c r is a positive integer with c\ — 1; that explains our earlier 
remark that A m (n) is a sum of more and more terms, all of degree less than q, as n increases. 
(Corollary 0, below, shows furthermore that min^ n r i = 0.) This justifies the description of 
A m (n) as given by finitely many polynomials of the same degree and leading coefficient over 
the whole range n > 0. 

One gets v m (n) immediately from X m (n) via Equation (pQ). 

To prove the sufficient condition for n to be large enough we employ the integral gain 
graph $. We can easily bound n ($). The maximum gain of any edge is the maximum 
abscissa magnitude of any (//, k) G Ap that could represent an attack on the m-row 
board, which means that \k\ < m. By symmetry it is sufficient to consider < k < m. If all 
qj < 1, no horizontal attack is possible so we can restrict k to be positive. As a path has no 
more than q — 1 edges, all path gains are < (q — 1) • max \fi\. □ 

The strip interpretation of no(<&) is that it is the greatest breadth that can be attained, 
within the available height of m rows, by a sequence of moves that hits each row j at most 
qj times. 

We can find no($) exactly for some kinds of piece, including queens, bishops, and nightrid- 
ers. We suppose the nonhorizontal moves are all contained in the upper and lower of the 
four quadrants formed by the lines y = fax, where a -1 is a positive integer b. That means 
any move (x, y) satisfies < b if y ^ 0. (P may have bounded or unbounded hori- 

zontal moves.) For instance, a = b^ 1 = 1 for a queen or bishop. We further suppose that 
(x, ±bx) G Ap for every nonzero integer x. Finally, we assume there is exactly one piece in 
each row, so q = m. 

Proposition 3. With P as described and with all qj = 1, v m {n) is a polynomial of degree 
m on the domain n > n (&) = b\_ m ~ 2 \ but not on n > n (&) — 1. 

Proof. The gain graph has vertex set {v±, . . . ,v m } and edges (±bh)ei t i + h for i,i + h G [m] 
and h ^ 0. We want the maximum gain of a path; we do that by characterizing the paths 
of largest gain. Let W be such a path. Note that V is ordered by subscript and that, by 
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maximality, every edge in W is of the form \bh\ei^ + h- Thus we need only consider such edges, 
omitting any others; we may then simplify the notation to ei i+h- 

There are three kinds of internal vertices of W . At a zig, both neighbors are higher vertices. 
At a zag, both neighbors are lower. At a smooth vertex, one neighbor is higher and the other 
is lower. There are also the two end vertices of W and there might be isolated vertices, not 
in W. 

If Vk is a smooth vertex, it is in edges e^, of W with i < k < j or i > k > j. Replacing 
£ik e kj by eij leaves the gain unchanged. Thus, we may assume W has no smooth vertices. 

Suppose that after eliminating smooth vertices there is a vertex £ V(W). We can 
extend W from an endpoint V{ by the edge with gain b\k — i\. That increases the gain, 
contradicting the maximality of W. Consequently, no Vk can exist, so = V, and 

indeed no path with maximum gain can have a smooth vertex. 

All zigs are lower than all zags. If not, say k < k', Vk is a zag, and vy a zig. Thus W 
contains paths e^e^- with i,j < k and e^/e^y with > k! . Replace these paths with 
e-ik'&k'j and ev^ky- This increases the gain. 

There are now two cases depending on the parity of m. The endpoints Vi, Vj of W are 
both zigs, or both are zags, if m is odd, and one is of each kind if m is even. 

Suppose m is even. Extend W by the edge e^. Now we have a cycle C in which every 
vertex is a zig or a zag. All zigs have lower indices, so they are Vi, . . . , f r _i, v r where r = m/2; 
and the zags are v m ,v m -x, . . . ,v m+r+ i. It is easy to calculate that <p(C) = b[A(V\ + m] — 
bm 2 /2. We recover W by removing the edge of least gain, which is e rjr+ i of gain 6; thus, 
n Q ($) = b{m 2 -2)/2. 

If m is odd, we may suppose v^,Vj G L. If there is Vk G L with i < k, then interchanging 
Vi with Vk in W increases the gain. It follows that the zigs are v i, . . . , V( OT -3)/2, vt = f ( m -i)/2, 
and Vj = f( m +i)/2- Therefore, 



n ($) = <p(W) = b 

which simplifies to b(m 2 — 3)/2. □ 

Suppose the inverted slope a -1 is not an integer. Then n ($) < a -1 L ™^" 1 ] because the 
gain of 6jj is the greatest width attainable with height j — i, that is, L a_1 (j ~ *)J • Thus, the 
gain of an edge in a path of maximum gain is sometimes < — i). (This unavoidable fact 

requires some proof, which we omit.) A consequence is that u m (n) may be a polynomial on 
a domain including n > a" 1 [ m 2 ~ 4 J — 1, contrary to Proposition |3j Nevertheless, we obtain 
a more general bound than that of Theorem [TJ It is no longer necessary to assume a piece 
appears in every row but we do assume any row has at most one piece; thus q, the number 
of pieces, is also the number of occupied rows. 

Proposition 4. If all nonhorizontal moves are contained between the lines y = ±ax where 
a > 0, and if all qj < 1, then v m (n) is a polynomial of degree q on the range n > or x [ m 2 ~ 2 j . 



/ro+JA / m-L \ >> 



The proof is similar enough to that of Proposition |3] that it can be omitted. 
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5. The method of deletion and contraction 

To calculate the integral chromatic function we need vertex weights and the formula for 
deletion and contraction of an edge in a weighted integral gain graph: 

(3) *?*,/») W = X%,h)\e(n) ~ X%, h )/e( n )- 

This formula is valid for any link e, by [21 Theorem 3.4], and translates into a formula for 
v m (n) by Equation JT]). 

Another useful formula is multiplicativity. If $ is the disjoint union of $ 1; . . . , $ r , then 

z _ z z 

A($,/i) — X{<S>i,h) ' ' ' X(<S> r ,h)- 

Now we begin to compute X%ny ^ n our calculations the notation x + means the positive 
part of the real number x, that is, 

x + := max(a;, 0). 

The effect of a loop is simple. If it has nonzero gain k, it can be discarded because it 
stands for a constraint Xj ^ xi + k, which is always satisfied. A loop with zero gain, however, 
forces xK 77 -) = f° r a U n because the corresponding constraint is unsatisfiable. 

Example 1 (A Vertex with Loops). If $ has one vertex Vi, all edges are loops. Thus, 

(n — hi) + if there is no zero loop, 



X(<3>,h) [ 



n 



if there is a zero loop. 



Example 2 (Isolated Vertices with Loops). If $ has no edges other than loops, then every 
component has a single vertex. By Example [1] and multiplicativity of the integral chromatic 
function, 

z , , I (n — • • • (n — if there are no zero loops, 

I q there is a zero loop. 

This gain graph is the basic type to which we reduce all other gain graphs by deletion and 
contraction. 

Example 3 (A Multiple Edge). Suppose $ has two vertices and all its edges have the form 
ei = ^iViV2, no two gains being equal. Let M = {//j} be the set of gains. If we contract 
one edge, the other edges become loops with nonzero gain, so may be ignored. The switched 
weights are (hi + fa, /12) if fa > and (hi, hi + \ fa\) if fa < 0; the contracted weight is the 
greater of the two switched weights. Thus, the deletion-contraction formula gives 

Xf*, h )(n) = (n-h 1 ) + (n-h 2 ) + 

(n - max(/ti + fa, h 2 )) + + ^ (n - max(/ii, /i 2 + f • 

i-.Hi>0 r.fii<0 ) 

If the original weights are both 0, this simplifies to 

x\ m {n) = (n+f-Y+n-\fa\y. 
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Figures [I] and [2] give examples of the deletion-contraction process symbolically, in the form 
of graph equations, and show the details of the contraction process with its two stages of 
switching (except for a zero edge) and contracting. Symbolically, Q 2 is the difference between 
the deletion and the contraction of e, and so forth. 



V 2 V 2 




v i v i 



Ql B 2 Q 2 /e 



Figure 1. The weighted integral gain graph Q 2} of order 2 with edges having 
gains —1, 0, 1 from vi to v 2 and vertex weights all 0, reduced by deletion and 
contraction of e = 0t>it> 2 . B 2 denotes Q 2 \e. 

Example 4 (A Small Problem of Queens and Bishops). To further explain the contraction- 
deletion method, we treat in detail the cases of queens and bishops with m = 2 and q± = 
q 2 — 1, that is, one piece in each of two rows, so that the (weighted integral) gain graph has 
order 2. Since there is one piece per row, u m (n) = x|( n )- 

The gain graph for each has vertices t>i and v 2 , one for each row, and all vertex weights 
equal to 0. The gain graph Q 2 for queens (see Figure [IJ has edges 0viV 2 , lviv 2 , —lv\V 2 
corresponding to the constraints x 2 ^ X\, x 2 ^ x\ + 1, x 2 ^ x\ — 1 because a queen going 
up one row moves 0, 1, or —1 squares to the right. The gain graph B 2 for bishops (Figure 
|2fa)) has the same edges but without the one with zero gain since a bishop cannot move 
vertically. From Example [3] we get the formulas 

(4) v B2 ( n ) = xl 2 {n)=n 2 - 2(n - 1)+ 

(5) v Q2 {n) = 4» = n 2 - n - 2(n - 1) + . 

In the contraction-deletion process, the first step is to reduce the queens graph to the 
bishops graph by deleting and contracting 0v\v 2 . Figure [T] shows this reduction. The loops, 
which do not influence the integral chromatic function, can be deleted. We deduce the 
formula 

X0»=xl>)-n, 
because the contracted weight is max(0,0). 

The next edge to go is lv\V 2 . In order to contract it we have to switch by r][ — 1, r)' 2 — 0, 
giving the switched graph B 2 shown in Figure Efb). Then we can contract, as shown in 
the figure. The weight in the contraction is hi 2 = max(l, 0) = 1 because the corresponding 
variable x% 2 has to satisfy both bounds x\ > 1 and x 2 > 0. We deduce that 

x! 2 (n)=X^ Vlt >)-(n-l) + . 

The notation —lviV 2 here signifies a graph with two vertices and one edge whose gain is — 1 
from v\ to v 2 , or +1 in the opposite direction. 
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Figure 2. (a) The weighted integral gain graph B 2 , of order 2 with edges 
having gains —1,1 from v\ to v 2 and zero weights, reduced by successive dele- 
tion and contraction of e' = lv\V 2 , then e" = —lv\V 2 . (b) B 2 switched by 
{V11V2) = (1, 0) so e' can be contracted, (c) The contraction process for e", 
switching B 2 \ e' by (r)'{,T)' 2 ) = (0, 1) so e" has gain 0, then contracting e". 



We delete and contract —\v\Vi in Figure Eta) and (c). Switching by rj" = 0,7/2 = 1 (since 
the positive direction of the edge is from t> 2 to V\) gives the gains and weights in Figure 
E^c), whence comes the symbolic reduction formula in the figure (again, the loops should be 
ignored) and thence the formula 

X-2 VlV2 (n)=n 2 -(n-l) + . 

Combining the steps we obtain @ , as we ought. 

The contraction-deletion formula implies a slightly more precise description of the inte- 
gral chromatic function, and thence of the counting function u m (n), than that stated in [2J 
Corollary 3.2], quoted earlier at Equation (T5]). The improvement comes from giving a new 
proof of based on deletion and contraction. 
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Corollary 5. The integral chromatic function xf$o)( n ) °f a nonempty O-weighted integral 
gain graph has the form of Equation fl2]) where the n ri are nonnegative integers with all 
n u = 0. 

Proof. In 03]), deletion preserves the order of the graph and the coefficient sign, while con- 
traction changes the sign and the parity of the order. Thus, after eliminating all links by 
repeated deletion and contraction, we have an alternating sum of integral chromatic func- 
tions x%r- h ) °f hnkless integral gain graphs. The sign of a term depends on the number of 
vertices lost in contraction. Only one of the graphs has q vertices; that is the one obtained 
by deleting every edge, and in it the weights remain as in the original gain graph. It follows 
that 

(6) xf m (n) = J2(- 1 ) q ^xf^ hj) (n), 

3 

summed over all the linkless gain graphs that result from the deletion-contraction process, 
where dj is the order of tyj. Every one of these graphs has at least one vertex so dj > 0. 

The way weights contract ensures that the weights in each tyj are nonnegative, and also, 
since the original weights are all zero, that the minimum weight in tyj is 0. The integral 
chromatic function of a linkless graph (ty, h) equals Y\ v&v {n — h v ) + unless it is identically zero 
(see Example [2]). The corollary follows by substituting this expression in <^ and collecting 
equal expressions. □ 

For sufficiently large integers n, the superscripts + can be dropped. When all the terms 
are multiplied out, there is no cancellation; the resulting formula is 

q k 

(7) xf$, )H = ^2(-l) q ~ d n d ^ c j<?d J -d(n jl , n jdj ) for n > maxn^, 

d=0 j=l 

where a^-d denotes the elementary symmetric function of degree dj — d (identically if 
d < or d > dj, 1 if d = dj), all Cj > 0, and all nji > 0. 



6. Generating functions 

As we have seen, the counting functions z/ m (n) that we get, and more generally those that 
arise in counting lattice points outside an integral affinographic hyperplane arrangement [2] , 
are integer-weighted sums of terms of the form 

(8) a n = (n - ni) + (n - n 2 ) + ■ ■ ■ (n - n r ) + 

where n\ < n 2 < ■ ■ ■ < n r . The generating function of such a term has a simple form, 
p(t)/(l-t) r+1 , where p(t) is a polynomial which we can calculate exactly. That enables us 
to write down the generating function of u m , which is 



N m (t) :=E^(n)t" = -^ 



(1 - f)9+! 
n=0 v ' 

where q is the number of vertices of the gain graph, i.e., the number of pieces on the board. 
This is a special case of the general behavior of the generating function of lattice-point 
counts of a weighted integral gain graph, X^ t f^(t) := J2^Lo xf* /t)( n X™' "which has the 
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form = P(ia,h) (£)/(! ~~ t) q+1 where the numerator is a polynomial. The deletion- 

contraction formula extends immediately to generating functions: 

for a link e; hence, the numerator polynomials satisfy 

P{%h)(t) = P(%h)\e(t) " (1 " t)P (%h)/e (t). 

Let us develop the generating function of a term of the form ([H]). Let di = n r — for 
1 < i < r and let be the kth elementary symmetric function of d±, . . . , d r _\. Then 

oo r . , \ 

(9) ^aX = f^(-ir, A,U) 



'r-J 



n=0 j=l v ' 

where Aj(t) is the Eulerian polynomial, i.e., the polynomial such that 



71=0 



(see [HE]). The first few Eulerian polynomials are A (t) = 1, Ai(t) = t, A 2 (t) = t 2 + t, 
A 3 (t) = t 3 + At 2 + t, and A A {t) = t A + lit 3 + lit 2 + t. 

The proof of is by substituting for n the variable p := n — n r and multiplying out the 
resulting product p(p + di) ■ ■ ■ (p + rf r -i), in which the coefficient of p> is s r _j. 

The most useful cases for our examples are the generating function of (n — ni) + , which is 



t 



ni+l 



(i-tr 

that of n(n — n 2 ) + , which is 

t" 2+1 [(n 2 + l)t- (w 2 - 1)] 

(l-f)» 

and that of n 2 (n — n 3 ) + , which is 

t" 3+1 [(n 3 + l) 2 t 2 - 2(n 2 - 2) + (n 3 - l) 2 ] 

(i-ty 



Example 5 (j4 Small Problem of Queens and Bishops, continued from Example \4§. The 
generating functions are, for bishops, 

A 2 (t) t 2 2t 3 -t 2 + t 

and for queens, 

2t 3 -t 2 + t t 2t 3 

N Q2 (t) - 



(1-tr il-t) 2 (i-ty 
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7. Some chess pieces, real and fairy 

Let us look at some more examples. We assume in all of them that m > 1 and there 
is one piece in each row (all qj = 1). We write v m (n) for the number of nonattacking 
arrangements with m rows, N m (t) := Yl^=o u m(ji)t n for its generating function, and v m (n) 
for the polynomial when n is large. In all the formulas, we assume n > 0. 

Every example has Ui(n) = n for n > 0, and Ni(t) —t/(l — t) m+1 . Also, all vertex weights 
are zero in the weighted integral gain graph for each piece. 

What we would most like would be to find a pattern in the generating function numerator 
for a fixed piece as the height m varies; but one would expect any regularity to show up 
only when the height is somewhat large because when it is small the moves are strongly 
constrained by the narrowness of the strip. The solutions we calculate are, unfortunately, 
too small to show a pattern. 

We calculated the function m several examples by reducing to gain graphs whose 

edges are all loops (and consequently may be discarded, unless a loop gain is zero, in which 
case the graph itself may be discarded), using deletion and contraction, as illustrated in 
Examples H] and [51 The relationship between the graph quantity xl an d I'm is that they 
are equal, according to Equation (pQ). We were able to do hand calculations for bishops and 
queens with m < 3 and knights with m < 4; we verified and extended these results by a 
computer package that carries out deletion and contraction of weighted integral gain graphs 
[I] with the appropriate simplifications; that is, it discards nonzero loops and graphs with 
zero loops. 

Example 6 (Rooks). Rooks, of course, are classically easy: v m (n) = (n) m , the falling factorial. 
The theorem says this is a polynomial for n > 0, which is true. The generating function is 
N m (t) = Y.n=o( n )mt n = t m m\/(l - t) m+1 , as is well known. 

Example 7 (Bishops). Proposition [3] applies with 6 = 1. Thus, v m {n) is a polynomial for 
n > i(m 2 — 3) but not if we extend the domain by one more integer, to n > |(m 2 — 5). 
The bishops gain graph B m has edges fiViVj with gains ±(j — i) for i < j. We found that 

u 2 (n) = n 2 - 2(n - 1)+, 

u 3 (n) = n 3 - n{A(n - 1)+ + 2(n - 2) + } 

+ {2(n - 1) + + A(n - 2)+ + 4(n - 3)+}, 
z/ 4 ( n ) = n 4 - n 2 {6(n - 1)+ + 4(n - 2)+ + 2(n - 3) + } 

+ n{4(n - 1)+ + 12(n - 2)+ + 16(n - 3)+ + 4(n - 4)+ + 4(n - 5) + } 

+ {4 [(n - 1) + ] 2 + 4(n - l) + (n - 3)+ + 4 [(n - 2) + ] 2 } 
- {2(n - 1)+ + 8(n - 2)+ + 34(n - 3)+ + 12(n - 4)+ 
+ 20(n - 5)+ + 4(n - 6)+ + 2(n - 7)+}. 
The polynomials for large n are 

i>2(n) = n 2 — 2n + 2 for n > 1, 

v^(n) = n 3 - 6n 2 + 18n - 22 for n > 3, 

u 4 (n) =n A - I2n 3 + 72n 2 - 234n + 338 for n > 7, 
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z/ 5 (n) = n 5 - 20n 4 + 200n 3 - 1192n 2 + 4132n - 6562 for n > 11, 

z7 6 (n) = n 6 - 30n 5 + 450n 4 - 4198n 3 + 25238n 2 - 91572n + 155220 for n > 17. 

The lower limits of validity can be read off from the piecewise-polynomial formulas (shown 
for m < 4; we computed m = 5,6 as well) and are as predicted by Proposition [3J The 

coefficients of n m ~ 1 and n m ~ 2 in the polynomials v m (n) are the values — 2r^j and 2(™) . 
The first can be explained intuitively by the fact that when n is large a bishop in row i 
blocks 2{m — i) squares in higher rows (see Proposition [6]). We have no explanation of the 
second. 

Generating functions are 

, . 2t 3 - t 2 + t .. 4t 6 - 4t 5 - 2t A - t 3 + Tit 2 - 13t 

N ^t) = —R Iv^-> ^ 3 (t) 



y l-t) 3 ' dV ' (l-^) 4 

2t n - 2t 10 + Ut 9 - 38t 8 + 54t 7 - 74t 6 + 64t 5 - 43t 4 + 213t 3 - 247t 2 + 81* 



(I-*) 5 

As a further illustration of the method of gain graphs we give some details of the solution 
for z/ 3 . The gain graph B 3 , shown in Figure [S^a), has vertices f i, v 2, 1*3 with v iV j-edges having 
gains ±(j — i) for i < j. We delete and contract edges in the arbitrarily chosen order If if2, 
— If if 2, If 2V3, — If 2^3, after which each component is a vertex or a multiple edge. The 
treatment of these multiple edges is as in Example [31 so is not shown here. 

Example 8 (Queens). As with bishops, Proposition [3] applies with 6 = 1. Thus, v m {n) is a 
polynomial for n > |(m 2 — 3) but not for n > |(m 2 — 5). 

The queens gain graph Q m has edges /if jf j with gains \i — 0, ±(j — z) for z < j; it is the 
bishops graph B m together with edges. The results are: 

u 2 {n) =n 2 - {n + 2(n - 1) + }, 

v 3 (n) = n 3 - n{3n + 4(n - 1)+ + 2(n - 2) + } 

+ {2n + 8(n - 1) + + 8(n - 2) + + 4(n - 3)+}, 
z/ 4 ( n ) = n 4 - n 2 {Qn + 6(n - 1) + + 4(n - 2)+ + 2(n - 3)+} 

+ n{lln + 30(n - 1)+ + 32(n - 2) + + 26(n - 3)+ + 4(n - 4) + + 4(n - 5) + } 

+ {4 [(n - 1) + ] 2 + 4(ra - l) + (n - 3) + + 4 [(71 - 2)+] 2 } 

- {6n + 40(n - 1) + + 72(n - 2) + + 94(n - 3) + + 36(n - 4) + 

+ 32(n - 5) + + 4(n - 6) + + 2(n - 7) + }. 

The polynomials for large n are 

V2(ri) = n 2 — 3n + 2 for n > 1, 

^3(^) = — 9n 2 + 30n — 36 for n > 3, 

p 4 (n) = n 4 - 18n 3 + 139n 2 - 534n + 840 for n > 7, 

1*5(71) = n 5 - 30n 4 + 407n 3 - 3098n 2 + 13104n - 24332 for n > 11, 

P 6 (n) = n 6 - 45n 5 + 943n 4 - 11755n 3 + 91480n 2 - 418390n + 870920 for n > 17. 

The ranges of validity of the polynomials agree with Proposition [3j The second term in the 
polynomial has coefficient — 3(™) , which as with bishops is the sum of the number of squares 
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-1, lA 




-1, l/< 



v 2f |-2, 2 - "of -2, 2_ {-3,-1,1 
by(b) -1" 



-1,1" 



(a) 



-1, In 



-2,2 _ 



by (c) 




m 

{-2,0,2 _ {-3,-1,1 
7 



-1a 



by (d) 



-2, 2 _ i,-2, 2 _ m-2,0,2_ h-3,-1, 1 



by (e) 



2, 2 _ A-1,3 _ / V -2, 2 _ m-2,0,2_ a— 3, —1, 1 



FIGURE 3. (a) The weighted integral gain graph for 3 bishops in 3 rows, 
showing the contraction process for (b) the first, (c) the second, (d) the third, 
and (e) fourth steps of deletion and contraction, e denotes If if 2, e' = —If if 2, 
and e" = lv 2 v 3 . For simplicity, parallel edges are drawn as a single edge with 
multiple gains. Nonzero loops are omitted. 

excluded by each piece in rows higher than itself. The third coefficient does not follow an 
easily identifiable pattern. 

The generating functions of the exact numbers are 

2t 3 a r / \ 4 * 6 - 4 * 4 - 4 ^ 3 + 24t 2 - Ut 

N 2(t) = 7T-^> N *(t) = 



^4(^) 



(1-O 3 ' (i-*) 4 

2t u - 2t w + 26t 9 - 50t 8 + 78t 7 - 130t 6 + 66t 5 - 42t 4 + 484t 3 - 640t 2 + 232t 



Example 9 (Knights). A widest knight's path lets the knight rise from the bottom row to 
the top, one row at a time, moving 2 steps to the right at each row. There are m — 1 moves 
that cover a width of 2m — 2. We conclude that v m (n) is a polynomial for n > 2m — 2 but 
not for n > 2m — 3. 
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B 3 - ee'e" 



The knights gain graph Kt m has edges ±2v{Vi + i and ±lviV i+2 . 

We calculated the first few cases, m < 3 by hand and m < 6 on the computer, via gain 
graphs. 

u 2 (n) = n 2 -2(n-2) + , 

u 3 ( n ) = n 3 - n{2(n - 1) + + 4(n - 2)+} + {6(n - 2) + + 4(n - 3) + + 2(n - 4)+}, 
z/ 4 (n) = n 4 - {4n 2 (n - 1) + + 6n 2 (n - 2) + } 

+ { 16n(ra - 2) + + 12n(n - 3) + + 4n(n - 4) + + 4[(n - 1) + ] 2 + 4[(n - 2)+] 2 } 
- {I2(n - 2) + + 24(n - 3) + + 18(n - 4) + + 4(n - 5) + + 2(n - 6)+}. 
The polynomials for large n are 

v 2 {n) = n 2 - In + 4 for n > 2, 
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z7 3 (n) = n 3 - 6n 2 + 22n - 32 for n > 4, 

z7 4 (n) = n 4 - 10n 3 + 56n 2 - 168n + 220 for n > 6, 

z7 5 (n) = n 5 - 14n 4 + 106n 3 - 478n 2 + 1248n - 1480 for n > 8, 

z7 6 (n) = n 6 - 18n 5 + 172n 4 - 1028n 3 + 3956n 2 - 9154n + 9852 for n > 10. 

The ranges of validity, read off from the piecewise-polynomial formulas, agree with the 
preceding calculation. The size of the second coefficient, as with bishops and queens, is 
ci = 4m — 6, which is the number of squares excluded in higher rows as well-separated 
knights are inserted from the bottom row up, that is, 4 squares for each knight except the 
one in the next-to-last row, which excludes 2 squares in the last row, and the last knight, 
which excludes no squares since there are no rows below it. (See Proposition [6j) The third 
coefficient appears to be a quadratic function of m, = 8m 2 — 22m + 16. (As with bishops, 
we have not proved this.) 

As for generating functions, we have 

2t 4 - 2t 3 + t 2 + t 2t 7 - 8i 4 + 3i 3 + 24t 2 - 15t 

N>2.{t) — -r^ , N 3 {t) 



N 4 (t) 



(1-t) 3 ' aw (l-^) 4 

2t 10 - 2t 9 + 12t 8 - 20t 7 - 6t 6 + 18t 5 - 25i 4 + 181t 3 - 199i 2 + 63t 



Example 10 (Nightriders) . This fairy chess piece has the move of a knight extended indef- 
initely in a straight line; that is, it moves along lines of slope ±2 ±1 . Therefore Proposi- 
tion [3] applies with 6 = 2, and v m {n) is a polynomial for n > 2\m 2 /2\ — 2 but not for 
n > 2[m 2 /2\ - 3. 

The gain graph Nr m for nightriders has edges ±2(j — i)v(Vj if j — i is nonzero, and edges 
±|(j — i)v{Vj if j — i is nonzero and even. The former correspond to moves of any number 
of rows in the directions ±(±2, 1) and the latter to moves of an even number of rows in the 
steeper directions ±(±1,2). 

Using the computer, we calculated the same data as in the previous examples for m < 6. 
The counting functions for m = 1, 2 are the same as for the knight. For 2 < m < 4, they are 

u 2 (n) =n 2 -2(n-2) + , 

u 3 (n) = n 3 - n{2{n - 1)+ ± 6(n - 2)+} ± {l0(n - 2) + ± 4(n - 3) + + 2(n - 4)+}, 
i/ 4 (n) = n 4 - {4n 2 (n - 1)+ + I2n 2 {n - 2)+} 

+ {40n(n - 2) + + 16n(n - 3) + + 2An(n - 4) + 

+ 4[(n - 1) + ] 2 + 8(n - l) + (n - 2) + + I2[(n - 2)+] 2 } 
- {58(n - 2)+ ± 56(n - 3)+ + 100 (n - 4) + + 2A(n - 5) + + 24(n - 6) + }. 
The polynomials for large n are 

z/ 2 (^) = n 2 — 2n + 4 for n > 2, 

z>3(n) = n 3 — 8n 2 + 34n — 56 for n > 4, 

p 4 (n) = n 4 - 16ra 3 + 132n 2 - 566n + 1016 for n > 6, 

z/ 5 ( n ) = - 28n 4 + 390n 3 - 3100n 2 + 13600n - 25676 for n > 8, 
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P 6 (n) = n 6 - 42n 5 + 876n 4 - 10974n 3 + 84720n 2 - 374678™ + 730408 for n > 10. 

The second coefficient, up to sign, is 2(^) +2( L(m+ 2 1)/2j ) +2( [m 2 /2j ), as Proposition implies. 
For nightriders, as with queens, the third coefficient Ci does not have a discernable pattern. 
The generating functions are 

2t 4 - 2t 3 + t 2 + t 6t 7 - 8t 6 + 8t 5 - 16t 4 + 5t 3 + 32i 2 - 21* 

N2\t) — J-, 77^ , N 3 (t) 



v l-t) 3 ' w (l-^) 4 

24t 10 - 48t 9 + 100t 8 - 196t 7 + 178t 6 - 114t 5 + 15t 4 + 144t 3 - 585t 2 + 209t 

8. The polynomial 



We observed in the examples that the polynomial v m {n) which gives ^ m (n) for wide boards, 
i.e., for large values of n, has second term — c\n m ~ l where C\ is the sum over all rows i of 
the number of squares attacked in higher rows by a piece in row i. There is a general result 
here. Suppose there are pieces in row i. The polynomial can be written in the form 

n q — c\n q ~ x + c 2 n q ~ 2 + • • • 



v m (n 



with q = q 1 H V q m . 

Proposition 6. Assume that (0,0) is a move. Then the magnitude of the second coefficient 
ofv m (n) is 



E 

i=l 



2 ) + Qi 



j=i+i 

where is the number of squares in row j that are attacked by a piece in row i. When every 
row has one piece, 

m 

ci = /,a>i, 

where ai is the number of squares in rows % + 1, . . . , m that are attacked by a piece in row i; 
thus, C\ equals half the maximum total number of squares attacked (nonhorizontally) by one 
piece placed in each of the m rows. 

One calculates the number of squares attacked ignoring any limitations due to inadequate 
board width. 

Proof. The proof is by inclusion and exclusion using labelled pieces. Let w be the greatest 
width of a nonhorizontal move that is possible on a board of height m, and assume n is very 
large compared to w. Number the pieces P%, P2, . . . , P q so that the pieces in higher rows have 
higher numbers. 

Imagine the board extended indefinitely to left and right by "imaginary" squares supple- 
menting the "real" squares of the m x n strip. Let U be the set of all placements of the q 
pieces, attacking or not, and let A^i be the set of placements in which P^ attacks Pi and the 
latter is in a higher row. Let i be the row of P^ and j that of Pi. Then 

\U\ = Y\_( n )qh = n q — ( + terms of lower degree, 

h=l h=l ^ ' 
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where (n) r is the falling factorial. 

To calculate \A^\ we place all pieces except Pi, and then we put Pi into any square (in 
row j) attacked by Pk, including occupied and imaginary squares. Call the set of all such 
arrangements, in which Pk attacks Pi and the other pieces may be attacking, A'. Some of 
these arrangements have Pi in an imaginary square; call the set of those arrangements B. 
Some of the arrangements have Pi in an occupied square; call the set of those arrangements 
C. Since B and C are disjoint, \A k i\ = \A'\ — \B\ — \C\. We estimate \A'\, \B\, and \C\ in 
order to estimate \A\. 

First, | A' | is the number of ways to place q^ pieces in each row h, except that only — 1 
pieces appear in row i, and then to put Pi in any of the real or imaginary squares in row 
j that are attacked by Thus, 

1^1 = i n )qi ■ ■ ■ ( n )qj-i( n ) q] -i( n )q ]+ i • ■ • ( n )q m a ij = aij^' 1 + lower terms. 

Second, \B\ counts only placements where Pk lies within distance w from the ends of the 
strip, because it must be able to attack an imaginary square. If we count such arrangements 
by placing Pk and Pi last, there are no more than n q ~ 2 ways to place the other pieces in real 
squares and then at most 2w ways to lay down Pk and ways to put down Pi. That is, \B\ 
is bounded by a polynomial in n of degree less than q — 1. 

Third, |C| counts placements where Pi is attacked by Pk and overlies another piece in its 
row. Let us again lay down Pk and Pi last. The number of places for Pk is not more than 
a>ij(qj — 1), because Pk has to be attacked by one of the qj — 1 pieces already placed in row 
j, and because aji = aij by symmetry of moves. Then the number of places for Pi is also no 
more than qj — 1. Thus, |C| is bounded by a polynomial in n of degree less than q — 1. 

It follows that \A k i\ = a^n q ~ x + a quantity bounded by a polynomial of degree at most 
n q - 2 . 

The size of the intersection of two sets of bad placements, \Aki PI Afc/j/|, is bounded by the 
number of placements in which Pi and Py are attacked by Pk and P&, respectively, which 
is in turn bounded by a^aj/ j>n q ~ 2 , a polynomial of degree less than q — 1. Likewise, the 
intersection of more than two sets Aki is bounded in size by a polynomial of degree less than 
9-1- 

The proposition now follows from the formula of inclusion and exclusion combined with 
the fact that z/ m (n) is known to be a polynomial of degree q. Taking only the terms of highest 
order, and letting (k, I) range over pairs of piece indices for which Pk is in row i and Pi is in 
row j with i < j, 



P| A ki \ = \U\ - \ Aki\ + lower terms 

k,l 

m f \ 

( In 9-1 — qiqjdijn q ~ 1 + lower terms. □ 



in 

,1 



h=l x / i<j<m 



The especial importance of C\ is shown by a probabilistic application^ 

Corollary 7. The 'probability that a random assignment of one piece in each row is nonat- 
tacking is asymptotic to 1 — C\n~ x . Assuming that (0,0) is a move, the probability that a 
random assignment of qt pieces in distinct positions in row i is nonattacking is asymptotic 
to l-n-^Y^iKjQiQjaij- 



1 We thank undergraduate student Christian Noack for raising this question. 
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Proof. The first probability is v m {n)/n m = 1 — C\n 1 + 0(n 2 ) 
The second probability is 

v m {n) 1 — ex??, -1 + 0{n~ 2 ) 



= [1 _ Cin -i + 0(n" 2 )] [l + n- 1 J2( \ + 0{n~ 2 ) 

i=i ^ ' 

by the geometric series (when n is sufficiently large), 

+ 0(n~ 2 ). □ 



1 -n~ l 



,2 

i=l 



We did not find a formula for C2. The solution for bishops and the proof of Proposition 
suggest there may be one but that it is not usually simple. 
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Note added February, 2010. After acceptance of this article, Vaclav Kotesovec (January 
25, 2010) informed us of the prior discovery of the polynomials z>fc(n) for queens (Example 
ED by Pauls [TD] (k = 3, 1874), Tarry [TJ] (k = 4, 1890) and Kotesovec [3 [8] (k = 5,6,7, 
1992). Initial values of the polynomials can be found in the On- Line Encyclopedia of Integer 
Sequences [TTJ, sequences A061989-A061993. The references were provided by Kotesovec. 
Kotesovec's Web site [9] has many other formulas for placing non-attacking chess pieces. 
Furthermore, Kotesovec has produced (February, 2010) a conjectural simple formula for the 
third coefficient of the queens polynomials. 
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